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Within the framework of quantization of the macroscopic 
electromagnetic field, equations of motion and an effective 
Hamiltonian for treating both the resonant dipole-dipole in- 
teraction between two-level atoms and the resonant atom-field 
interaction are derived, which can suitably be used for study- 
ing the influence of arbitrary dispersing and absorbing mate- 
rial surroundings on these interactions. The theory is applied 
to the study of the transient behavior of two atoms that ini- 
tially share a single excitation, with special emphasis on the 
role of the two competing processes of virtual and real photon 
exchange in the energy transfer between the atoms. In par- 
ticular, it is shown that for weak atom-field interaction there 
is a time window, where the energy transfer follows a rate 
regime of the type obtained by ordinary second-order pertur- 
bation theory. Finally, the resonant dipole-dipole interaction 
is shown to give rise to a doublet spectrum of the emitted 
light for weak atom-field interaction and a triplet spectrum 
for strong atom-field interaction. 

PACS numbers: 42.50. Ct, 42.50.Fx 42.60.Da, 80.20.Rp 



I. INTRODUCTION 

Recently, several implementations of quantum logic 
gates relying on the resonant dipole-dipole interaction 
between two atomic qubits have been proposed. The 
atomic qubits could be impurity atoms in the condensed 
phase M or quantum dots embedded in a semiconductor 
Here and for many other applications in practice, the 
question of the influence of real material surroundings on 
the mutual interaction of the dipoles arises. Moreover, 
tailor-made material surroundings offer the possibility 
of controlling the mutual interaction of the dipoles. In 
particular, non-absorbing photonic crystals || and non- 
absorbing Fabry-Perot cavities [Q,[5| have been studied 
within the framework of mode expansion. Effects of ma- 
terial dispersion and absorption have been taken into ac- 
count for bulk material , Fabry-Perot cavities B , and 
microspheres 0. 

In the regime of weak atom-field coupling, the mutual 
resonant interaction of atoms has typically been char- 
acterized by an effective two-body potential involving 
atomic variables only. By analyzing a onc-dimcnsional 
cavity model, it has been shown that this concept may 
fail to give a correct description of the interaction at least 
in the strong-coupling regime, where the electromagnetic 
field degrees of freedom can no longer be eliminated M. 



In that case, it may be instructive to treat the interac- 
tion of the atoms with the on-resonant part of the elec- 
tromagnetic field exactly and the interaction with the 
off-resonant part in a perturbative manner. In particu- 
lar, for two two-level atoms that resonantly interact with 
cavity-type modes, an effective Hamiltonian of the form 

H cS = h^uxai&x + Ti (u>u a o-a^A + u UB a B a B 

A 
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has been proposed [Gil. Here, the first term is the energy 
of the cavity modes[of frequencies ui\ and photon de- 
struction (creation) operators a\ (ajj], the second term 
is the energy of the two atoms, where &a(b) an d "a(b\ 
are the Pauli operators of the two-level atom A(B), and 
Tilj Ua{b) is the upper-state energy, the third term is the 
energy of the so-called resonant dipole-dipole interaction 
of the atoms to each other (TiM - coupling strength), 
and the fourth term is the familiar interaction energy 
between the atoms and the cavity modes in the rotating- 
wave approximation. Apart from the fact that the cou- 
pling parameters in the Hamiltonian (|l|) are not speci- 
fied and their relation to each other thus remains unclear 
(also see fll||ll"|j ), material absorption cannot be taken 
into account, because of the underlying concept of mode 
decomposition. 

In this paper we give a rigorous derivation of an effec- 
tive Hamiltonian, which - although at first glance looks 
like that in Eq. (Q) - applies to atoms in arbitrary dis- 
persing and absorbing material surroundings, with the 
coupling parameters being well defined. For this pur- 
pose, we start from the multipolar Hamiltonian govern- 
ing the motion of the coupled system of two-level atoms 
and the medium-assisted electromagnetic field, with the 
medium being described in terms of a spatially vary- 
ing permittivity that is a complex function of frequency 
(for a review, see 13 ). The quantity that essentially 



governs the strength of the atom-field interaction is the 
(classical) Green tensor of the inhomogeneous Helmholtz 
equation with the space- and frequency-dependent com- 
plex permittivity of the material surroundings. Treat- 
ing the off-resonant part of the atom-field interaction in 
a (coarse-grained) Markov approximation, whereas the 
on-resonant part of the atom-field interaction is left in 
the original form, we derive the density-matrix equations 
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of motion for the coupled atom-field system, where the 
off-resonant atom-field interaction is eliminated in favor 
of the resonant dipole-dipole interaction. We then show 
that under certain conditions an effective Hamiltonian 
that governs the motion can be constructed. The results 
show that the resonant dipole-dipole interaction (via the 
off-resonant atom-field interaction) and the on-resonant 
atom-field interaction, respectively, are essentially deter- 
mined by the real and the imaginary part of the medium- 
assisted Green tensor. Whereas the spontaneous decay 
is only determined by the imaginary part, the mutual in- 
teraction of the atoms is determined by both the real and 
the imaginary part of the Green tensor. 

To illustrate the theory, we examine, for both weak 
and strong atom-field interaction, the temporal evolution 
of two two-level atoms that initially share a single exci- 
tation, with special emphasis on the interatomic energy 
transfer. In particular for weak atom-field interaction, 
from the exact time dependence of the upper-state pop- 
ulation of the acceptor atom we infer an energy transfer 
rate. We show that it is essentially the same rate as the 
commonly used rate, i.e., the transition probability per 
unit time which is obtained by means of Fermi's golden 
rule in second-order perturbation theory with regard to 
the original Hamiltonian. Finally we address the problem 
of the influence of the resonant dipole-dipole interaction 
on the spectrum of the light emitted by the atoms in the 
limits of weak and strong atom-field interaction. 

The paper is organized as follows. In Section [0] the 
density-matrix equations of motion for the system that 
consists of the atoms and the resonant part of the elec- 
tromagnetic field are derived. The problem of deriving 
them from an effective Hamiltonian is studied in Section 
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The temporal evolution of two atoms that initially 
share a single excitation is studied in Section [Fv|. The 
problem of determining an energy transfer rate is con- 
sidered in Section and Section |v| is devoted to the 
spectral properties of the emitted light. Finally, a sum- 
mary and some concluding remarks are given in Section 
VII. 
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Here, f (r, uj) and V (r, io) are bosonic fields which play the 
role of the fundamental variables of the electromagnetic 
field and the medium, including a reservoir necessarily as- 
sociated with the losses in the medium, G(r, r', uj) is the 
classical Green tensor, and e(r, uj) — £R,(r, uj) + ie\(r, uj) 
is the complex (Kramers-Kronig consistent) permittivity. 
It should be pointed out that there are no direct Coulomb 
forces between particles in the Hamiltonian (|2|) ; all inter- 
actions are mediated by the medium-assisted electromag- 
netic field. 

With regard to the interaction of the atoms with the 
electromagnetic field, it is convenient to decompose the 
latter into an on-resonant part (denoted by J Q °° duj . . .) 

and an off-resonant part (denoted byf Q °° duj . . .). Let 
us now consider the temporal evolution of the system 
that consists of the atoms and the on-resonant part of 
the electromagnetic field. If O is any system operator, 
we may write its equation of motion in the Heisenberg 
picture as, on recalling the Hamiltonian (H), 



O = 
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where 



H s = fd 3 r[ dujhuj&(r,uj)f(r,uj) +J2fa A a A 
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II. DENSITY-MATRIX EQUATION 
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dui 



d A E(r A ,w)+H.c. 



(6) 



Let us consider N two-level atoms [positions r A , transi- 
tion frequencies v A , dipole moments (A=l, 2, —,N)] 
that resonantly interact with the electromagnetic field 
via electric-dipole transitions in the presence of dispers- 
ing and absorbing bodies. The corresponding multipolar- 
coupling Hamiltonian reads as p^ , p^ | 



H = d'r dujhujr{r,uj){(r,uj) 
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where 



Note that in Eq. (|5|) normal ordering is adopted such 

that ^(r^w) is on the right-hand side and E^r^cj) is 
on the left-hand side in operator products. According to 

Eq. (||), the operators E(r A , uj) and E {r A ,uj) in Eq. (||) 
are thought to be expressed in terms of the basic- variable 
operators f(r,w) and F(r, uj) respectively. It is not diffi- 
cult to see that f(r,w) obeys the Heisenberg equation of 
motion 

f(r, uj) = —iujf(r,uj) 



(7) 



2 



We now solve Eq. ([?]) formally, insert the result into 
Eq. (||), apply a (coarse-grained) Markov approximation 
to the slowly varying atomic variables in the time inte- 
grals in the off-resonant frequency integrals, and assume 
that the off-resonant free-field is initially (t = 0) prepared 
in the vacuum state. After some algebra, we arrive at the 
following equation of motion for the expectation value of 
the system operator O (Appendix |X|) : 



(O) 
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where the notation YI'a A' indicates that A^= A' . In par- 
ticular, when O is identified with the atomic operator a a, 
Eq. (|) yields 
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In Eqs. @ and (^), the system Hamiltonian H$ is de- 
fined according to Eq. (Q), with lu a being replaced by & a, 
where 



&A = Ll)A <5a*a, 



5a* a = S A , A - 5\, A , 



(10) 



(11) 



8^=^r^ dAimG ^ A ^ A \ (12) 
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and the resonant interatomic coupling parameters are de- 
fined according to 



A'\a^A> — & A' A' + ^A*A'> 



(13) 
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[V - principal value]. The notation A* (A 1 * ) means that 
d.A (dA') in Eqs. ( |l2] ) and (|T4h has to be replaced with 
its complex conjugate (d* A ,). Applying the Kramers- 
Kronig relation to the Green tensor, from Eq. ( fbi] ) we 
derive that 



J AA' 



?1£qC 2 

and Eq. (|l^) thus becomes 



'— d A Re G(r A , r A > , &a> ) &a> 



J AA'-> 
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d A ReG(rA,rA',0JA')d A >- (16) 



Accordingly, from Eqs. (tL2j) and (llll) it follows that 
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In dealing with problems of atoms embedded in media, 
the atoms should be assumed to be localized in some 
small free-space regions, so that the Green tensor at the 
positions of the atoms can always be written as a sum of 
the vacuum Green tensor G v and the reflection Green 
tensor G R . Due to the singularity of TLeG v at equal 
space points, Eq. ( |l7| ) actually applies to the reflection 
part only. The vacuum part, calculated in many text- 
books, can be thought of as being already included in 

UJ A . 

It should be pointed out that the single-atom frequency 
shift 8 a* a [Eq. (^lj) or Eq. (|l7|)] (which is a real quan- 
tity because of the reciprocity property of the Green ten- 
sor) differs from the frequency shift 8 A , A obtained in the 
rotating- wave approximation (see, e.g., Ref. in the 
term S A , A , which results from the counter-rotating con- 
tributions to the Hamiltonian. Note that when the atoms 
are in free space, then Eqs. ([ill) an d (13) reduce to the 
result derived in Refs. pi 
right-hand side in Eq. 




cally the second term 5 A , A 



Since the first term on the 
171) can also be obtained classi- 
is sometimes termed the 
quantum correction. It is interesting to note that this 
quantum correction just corrects the rotating-wave re- 
sult. For that part of the frequency shift which is caused 
by the presence of the macroscopic bodies (mathemat- 
ically, by the reflection part of the Green tensor), the 
quantum correction may safely be neglected. 

From Eq. (16) it is seen that when the transition fre- 
quencies Cj a and Q A i of two atoms A and A' are dif- 
ferent from each other, then the strengths of the reso- 
nant dipole-dipole coupling I^a-A'I and |#a'*a| are not 
symmetric with respect to A and A' . Only if the differ- 
ences \oj a — uj A > | are small compared with the frequency 
scale of variation of the Green tensor, this asymmetry can 
be disregarded and oj a and uj A i may be replaced by an 
appropriately chosen mid- frequency, say (Q A + ui A >)/2. 
Whereas for atoms in free space such an approximation is 
unproblematic even for relatively large frequency differ- 
ences, the situation can drastically change if the presence 
of macroscopic bodies gives rise to a highly peaked Green 
tensor that rapidly varies with frequency. 
Recalling the relationship 



(6(t)) = Tr[p(0)O(t)] 

= Tr[/3(t)6(0)] = Tr[g(i)6(0)], 



(18) 



where O is an arbitrary system operator, p is the density 
operator of the overall system, and g is the (reduced) 
density operator of the system, and making use of the 
cyclic properties of the trace, from Eq. (||) we derive the 
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following equation of motion for the system density op- 
erator in the Schrodinger picture (Appendix |X|) : 



\ i E { S A'A> (&A 
I A, A' 



<J' A GA'Q- OA'QO A 



<j A o- A >Q- 



cr A , QUA 



H.c. 



(19) 



In this equation both the resonant interatomic interac- 
tion and the resonant atom-field interaction are taken 
into account, without any restriction to the strength of 
the latter one. It is worth noting that this type of equa- 
tion cannot be derived from an effective Hamiltonian in 
general. 

In particular in the case when the resonant atom-field 
interaction is sufficiently weak, it can also be treated in 
a Markov approximation. Let g be the reduced density 
operator of the atomic system and let us assume that the 
on-resonant part of the electromagnetic field is initially 
also prepared in the vacuum state (i.e., there is no exter- 
nal driving field). By tracing out both the on- and the 
off- resonant medium-assisted field variables, we arrive 
at the following master equation for the reduced density 
operator of the atomic system (see Appendix pf): 
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where 
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If the differences between the atomic transition frequen- 
cies are small in comparison with the frequency scale of 
variation of the Green tensor, so that the relations 



S ± ~ fi* 

°A*A> — At A* 



J A'A* 



and 



^A'A< — ^A'A' 

are valid, then Eq. (^) reduces to 



(22) 
(23) 



Q = -^i^2^A [<TAz, Q] + i E 5a ' a ' 



A, A' 



o\oA' , Q 



\ E T A'A> (o\a A >Q~ 2a A 'Qar A + ga A a A '^j , (24) 



A, A' 



which is of the same form as the master equation obtained 
on the basis of Kubo's formula for the field correlation 
functions 0. 



III. EFFECTIVE HAMILTONIAN 

Let us return to Eqs. (|J) and (|l^) and restrict our 
attention to the case when the difference between the 
atomic transition frequencies is small in comparison with 
the frequency scale of variation of the Green tensor, so 
that Eq. @ holds. Then Eqs. (§ and © reduce to 



<°> = 4([<M" 



•A'Vj^OA 1 

A, A' 



and 



(Hs - E MA*A><r A <rA'j , Q 



A, A' 



(25) 



(26) 



respectively. Recalling the definition of H$, we see that 
the motion of the system is governed by the effective 
Hamiltonian 

//>'oo 
d 3 r J duifiuj^{r,uj)i(r,uj) 

+ ^Z\^A^Az - E h$A*A>VA°A' 



A.A> 



E 



dw 



d A E(r A ,w) + H.c. 



(27) 



with E(ta,(jj) being given by Eq. (Q). 

From the above it is clear that a Hamiltonian of the 
type ([!]) makes only sense if the atomic transition fre- 
quencies are sufficiently near to each other. In that 
case, the Hamiltonian ( p7j ) is the desired extension of the 
Hamiltonian (|l|). The Hamiltonian ( p7|) is remarkable in 
several respects. Firstly, it applies to atoms surrounded 
by arbitrarily configured, dispersive and absorptive me- 
dia. Secondly, it goes beyond the rotating-wave approx- 
imation. Thirdly, it contains a resonant dipole-dipole 
coupling energy that is explicitly expressed in terms of 
the medium-assisted Green tensor, according to Eq. (p^). 

The first term in Eq. (27) describes, as before, the free 
medium- assisted electromagnetic field energy. The sec- 
ond term is the energy of the free atoms, which takes 
account of the medium-induced single-atom transition 
frequency shift. As already mentioned, the third term 
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represents the energy of the resonant dipole-dipole inter- 
action of the atoms. Note that both the single-atom tran- 
sition frequency shift and the interatom resonant dipole- 
dipole interaction result from the off-resonant atom-field 
coupling. The energy of the resonant interaction of the 
atoms with the electromagnetic field is given by the 
fourth term, which is of course not only responsible, e.g., 
for the single-atom spontaneous decay in the regime of 
weak atom-field coupling and the Rabi-type oscillations 
in the strong-coupling regime, but it also determines, 
together with the resonant dipole-dipole coupling, the 
mutual interaction of the atoms. In particular, in the 
limit of a <5-like field excitation the atoms are resonantly 
coupled to, the fourth term in Eq. (27) reduces to the 
interaction energy in the Tavis-Cummings model |p2| . 
Note that within the Tavis-Cummings model the reso- 
nant dipole-dipole coupling between the atoms cannot 
be described, because this model does not take account 
of the off-resonant atom-field interaction. In practice the 
excitation spectrum of any real electromagnetic field is 
always continuous, and the inclusion in the Hamiltonian 
of the off-resonant atom-field interaction is thus quite 
crucial, because it may lead to observable effects. 

It is worth noting that, as can be seen from Eq. (16), 



the strength of the resonant dipole-dipole interaction is 
determined by the real part of the Green tensor at dif- 
ferent space points. Thus, it is affected by a surrounding 
medium in a quite different way as the spectral excita- 
tion density of the medium-assisted electromagnetic field, 
which is determined by the imaginary part of the Green 
tensor at equal space points. In particular, it may happen 
that a high (low) excitation density and thus an enhanced 
(reduced) single-atom decay rate is accompanied by a re- 
duced (enhanced) strength of the resonant dipole-dipole 
coupling p[ . 

Clearly, strong resonant dipole-dipole interaction can 
only be expected if the atoms are sufficiently close to each 
other. To give a simple example of the effect of material 
absorption, let us assume that the atoms are embedded 
in bulk material of complex permittivity s(uj) . Using the 
bulk-material Green tensor (see, e.g., jl3|]), from Eq. ([l6| ) 
we find in the short-distance limit 



5 A* A' = 
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d*,R d^R 
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R 



d* A d A 



(28) 



(R=i"a — va 1 )- In free space, Eq. (|28j) reduces to the well- 
known result that the resonant dipole-dipole coupling 
simply corresponds to the (near-field) Coulomb-type in- 
teraction. It is seen that the characteristic R~ 3 distance 
dependence observed in free space is not changed by the 
medium. In the long-distance limit, we find that 



8 A' A' = 



AirheoRuj 2 ,, 



d* A d A ' 



d* A H d A 'R 
~ R R~ 



X cos 



nn(QA>) 



Cja'R 



exp 



-ni(wA') 



tJA'R 



(29) 



[n(uj) =71r(u;) +in\{uj) = y/sfuJ)], i.e., the (harmonically 
modulated) R^ 1 dependence observed in free space 
is changed to an exponential decrease according to 
exp[— ni(£jA')&A'R/c] due to material absorption. Thus, 
material absorption can drastically reduce the resonant 
dipole-dipole coupling strength with increasing mutual 
distance of the atoms. 



IV. TEMPORAL EVOLUTION OF A TWO-ATOM 
SYSTEM 



Let us use the effective Hamiltonian (|27|) to study a 
system of two-level atoms (resonantly) coupled to the 
medium- assisted electromagnetic field and assume that 
initially the atoms share a single excitation while the field 
is in the vacuum state. In the Schrodinger picture we may 
write, on omitting off-resonant terms, the state vector of 
the system in the form of 

|^(t))=^C A (t)e-^--)*|^)|{0}) 



d J r / duC Li (r,u,t)e 
o 



\L)f J (r,co)\{0}) (30) 



(u)= |2a^)- Here, \U a) is the atomic state with the 
Ath atom in the upper state and all the other atoms in 
the lower state, and \L) is the atomic state with all atoms 
in the lower state. Accordingly, |{0}) is the vacuum state 
of the rest of the system, and fj(r,u>)\{0}) is the state, 
where a single quantum is excited. 

From the Hamiltonian (p7|), the equations of motion 
for the slowly varying probability amplitudes Ca and C l 
read as 



C A (t) = iSA.A'e***-**^*® 



A' 
A'^A 

V^eon. Jo c 



oo , ,2 



x d* A G{r A , r, w) C L (r, w, t)e 



(31) 



1 ^ n — V 



£ d A > G* (t a , , r, w) Ca< (t)e^-^> )* . ( 32 ) 



By formally integrating Eq. (|32|) under the initial con- 
dition that Cz,(t = 0) = 0, and substituting the formal 
solution into Eq. (^lj), we obtain the following system of 
integrodifferential equations for the Ca- 



C A (t)= iSA-A'e^-^CA'it) 



A' 
A'^A 
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where 

K AA ,(t,t';u) 
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__g— uiA)t e i{u— &A')t' 

c 2 



hneo 

xd A lmG(r A ,r A ,,uj)d A 



(34) 



It should be pointed out that equations of the type (|3(]| ) 
- ( p4] ) can also be derived on the basis of the origi- 
nal Hamiltonian @ in the rotating wave approximation 

(Jq 00 du> . . . — > dui . . .), without the restrictive condi- 
tion that the Green tensor does not change on a frequency 
scale defined by the differences of the atomic transition 
frequencies Clearly in such an approach, the con- 
tributions of the counter-rotating terms to the frequency 
shifts and interatomic coupling strengths are disregarded. 
In order to get insight into the atomic motion on the basis 
of closed solutions of Eq. ( |33| ) , let us consider two atoms 
and restrict our attention to the limiting cases of weak 
and strong atom-field coupling. 



A. Weak atom-field coupling 

In the weak coupling regime, the integral expression 
in Eq. (^) can be treated in a (coarse-grained) Markov 
approximation, i.e., C A >(t') is replaced by C A i(t), and 



the time integral J Q dt' i 



,)(t-t') 



is replaced by the 



^-function (,(u) A i— uS). Since the frequency integral only 
runs over the resonance region, the (-function effectively 
acts as the ^-function, and we arrive at the following 
equations for the (slowly-varying) upper-state probabil- 
ity amplitudes of two atoms A and B: 

C A (t) = -\T A , A C A {t) +K, A * B e^ A - Q ^ t C B {t), (35) 

C B (t) = -\T B , B C B (t) +IC B , A e- i( - QA - QB ^C A {t), (36) 



where 



hF AB + i5 A i 



d A G(r A ,r B ,LD B )d B 



(37) 



We now make the simplifying assumption that the 
transition frequencies of the two atoms are nearly equal 
to each other, ub A ~ uj b , but allow for T A * A ^T B * B . The 
latter may happen, e.g., when the atoms A and B have 
different dipole matrix elements and/or different dipole 
orientations. It is then not difficult to solve Eqs. ([55]) and 
( |36| ) analytically. In particular, if the atom B is initially 
in the lower state, C B (t — 0) = 0, we obtain 



[±(T A * A -T B * B )+D] e D -^} 



C B = (>+*/2 - e - 1 ' 2 

where 



(38) 
(39) 



D = [\(T A * A - T B , B ) 2 + UC A , B K, B , A ] 1/2 , (40) 



D± = -UT A * A +T B . B )±D. 



(41) 



When the two atoms are identical and have equiva- 
lent positions and dipole orientations with respect to the 
material surroundings such that the relations 



F A > A — Ts'S, 
T A * B — T B * A , S A * B = 5 B * A 



(42) 
(43) 



are valid, we have K A * B —IC B * A , and T A * B , T B * A , 6 A * B , 
and S B * A are real quantities due to the reciprocity of the 
Green tensor. Then Eqs. (^) and ([39]) yield the following 
expressions for the upper-state occupation probabilities 
P M B){t)=\C A(B) (t)\ 2 : 



PA(B){t) = \ [cosh (IV B t) + (-)cos(25 A »Bi)]e- 



(44) 



It is worth noting that Eqs. ( p8| ) - (|4l]) [and thus 
Eq. (Q)] are valid for arbitrary dispersing and absorb- 
ing material surroundings of the atoms. In particular, 
substituting in Eqs. ( [H| ) and ( [37|) for the Green tensor 
the vacuum Green tensor, Eq. (|44|) reduces to that one 
obtained in Ref. jL6j. Accordingly, using the Green ten- 
sor for absorbing bulk material, the result in Ref. Q is 
recognized. 

From Eq. (|44|) a damped oscillatory excitation ex- 
change between the two atoms is seen. For sufficiently 
small times, Y B * B t-^il, and strong resonant dipole- 
dipole coupling, \5 A ' B \~^> T B * B , the oscillatory behav- 
ior dominates. This can typically be observed when 
the atoms are sufficiently near to each other, but can 
also be realized for more moderate distances with the 
interatom coupling being mediated by high-Q medium- 
assisted field resonances ||. In the opposite limit of 
weak resonant dipole-dipole coupling, P A {t) decreases 
monotonously while P B (t) features one peak which sepa- 
rates the regime of energy transfer from atom A to atom 
B at early times and the subsequent decay of the excited 
state of atom B. 



B. Strong atom-field coupling 

For the sake of transparency, we again consider identi- 
cal atoms that have equivalent positions and dipole ori- 
entations with respect to the material surroundings so 



G 



that Eqs. (42) and ( |43| ) hold. Introducing the probabil- 
ity amplitudes 

C±(t) - 2-5 [C A (i) ± C B (t)] e =F«A.Bt (45) 

of the superposition states 

|±)=2- 1 /2 ( | C / a)± | [ / b)); (46) 



where 



from Eqs. ( p3[ ) we find that the equations for C + (t) and 
C-(t) decouple, 



t p'oo 



C±(t)= At' duK ± (t,t';u;)e^* B ^- t ^C ± {t'), 
Jo Jo 

(47) 

where 

K ± (t,t';w) = K A * A (t,t';w)±K A , B (t,t';u). (48) 

Let us restrict our attention to the case when the ab- 
solute value of the two-atom term K A * B (t, t'; w) is of the 
same order of magnitude as the absolute value of the 
single-atom term K A * A (t,t';ui), so that there is a strong 
contrast in the magnitude of K + and K- . Typically true 
when the atoms are close to each other, this may also 
take place at interatomic distances much larger than the 
wavelength, e.g., as in the case of atoms situated at dia- 
metrically opposite positions near a microsphere ||. As 
a consequence, the strong-coupling regime is applicable 
to either the state |+) or the state |— ), but not to both 
at the same time. Assuming that the field resonance 
strongly coupled to the atoms has a Lorentzian shape, 
with oj m and Aw m being the central frequency and the 
half width at half maximum respectively, we can perform 
the frequency integral in Eq. ( |47j ) in a closed form, on ex- 
tending it to ±oo, 



/' 

Jo 



r± 



dujK±(t,t';ui) ~ -— (Aiv m ) 2 e 



\2 -i{u> m -ti A ){t-t' ) 



27T 



dtu ■ 



-i(uj—LU rn )(t—t ) 



r± 

2 



(CJ - LU m ) 2 + (Auj m ) 2 



(49) 



where 



A* A 



±r 



A*B, 



(50) 

m: 

with ui m in place of Cj a (~uj b ). Substituting Eq. (4£) 

(J47|), and differentiating both sides of 1 
suiting equation with regard to time, we arrive at 

C±(t) + [i (u m - Cj a ± (5 A * b) + Aw m ] C±(t) 



with T A * A and r ' a *b being defined according to Eq. (21) 

in p' ' " ' ^ 

into Eq. (|47|) , and differentiating both sides of the re- 



+ (fi±/2) 2 C±(i) = 0, (51) 



n± = v /2r±Aw r , 



(52) 



In particular for exact resonance, i.e., u> m ~ Cj a T $A* b , 
we derive 



C±{t) = 2-2 e -A^ m t/2 C os(fi±i/2) 



(53) 



(f2± S> Aw m ). For the probability amplitudes of the re- 
maining states |=f), which are weakly coupled to the field, 
we obtain 



c T (t) = 2-h~ r ^ 2 

(f7zp <C Auj m ). It then follows that 



(54) 



P 



A(B) 



(*) = ! 



,-r T t 



,-Au m t 



cos 2 (0±t/2) 



+ (-) 2e- (A ^ +r ^ t/2 cos(n±t/2) cos(25 A « B t)\ . (55) 

Here the upper (lower) signs refer to the case where the 
state |+) (|— )) is strongly coupled to the medium-assisted 
field. 

Even if damping is ignored, P A (b) (t) is not strictly 
periodic in general, because fl± and 8 A * B are not nec- 
essarily commensurate with each other. Roughly speak- 
ing, between three cases of approximately periodic mo- 
tion may be distinguished. 

(i) 4|5 A . B |»0±, i<2/fi ±: 

P A (t) = cos 2 (5 A * B t), 
P B {t) = sin 2 {5 A , B t). 

(ii) 4|^. fl |~n ±) t«i/|2a A . B -afc/2|: 

P A (t) = i [l + 3cos 2 (f>±t/2)] , 
P B (t) = \ [l-cos 2 (f! ± t/2)] . 

(iii) 4|5 A , B |<n±, i«l/|2* A . B |: 

P A (t) = cos 4 (ft ± i/4) , 
P B {t) = sin 4 (fi±t/4). 



(56) 
(57) 



(58) 
(59) 



(60) 
(61) 



From Eqs. ([36|) - (|6lj) the following time-averaged proba- 
bilities P A , P B , and Pj, = 1 — P A — P B are obtained, with 
the respective time integral being taken over one cycle, 
(i) P A ^P B ~ h and Pl=0. The excitation energy is 
periodically exchanged between the two atoms through 
virtual field excitations exclusively, (ii) P A = | , P B = | , 
and Pl = § ■ The two exchange channels - one channel 
through virtual and the other one through real field ex- 
citations - compete with each other and destructively 
interfere, leading to a partial trapping of the excitation 
energy in atom A. Note that this kind of energy transfer 
suppression cannot be observed, if the interference effect 
is disregarded as in Ref. Q. (iii) P A = P B = §, Pl = §• 
This case is typically observed in the long-distance limit, 
when the interatom energy exchange is dominantly me- 
diated by real field excitations. 
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V. RATE REGIME 

Let us return to the case of weak atom-field coupling. 
If the resonant dipole-dipole interaction is also weak, 
then, as already mentioned in Subsection IV A, the en- 



ergy transfer is one-way; that is, from atom A to atom B. 
Clearly, the efficiency of such a transfer regime is low, be- 
cause only a small portion of energy is passed on to atom 
B. In this case, a transfer rate w\ can be introduced 
according to 



!"1 



dPij(t) 



dt 



(62) 



where to is determined from the conditions that 
d 2 P B (t) _ n dP B (t) 



dH 



= 0, 



Co 



dt 



> 0. 



Note that this rate is not much different from that defined 
as the ratio between the maximum value of Pg (t) and the 
time belonging to it From Eq. (|3|) [together with 

Eqs. and ( f4l| ) for negligibly small Ka*b and K b *a 
thereinTit then follows that 



\K B < 



D 2 
1 , f 1 



A\ c D-tg 



D+e 2Dt _ p + +£) _) e Dto] 



(64) 



D V AD 



'+ 



2D[(D + + D_) 2 +4D + D_] 1/2 }^j, (65) 



where 



d = ±|r A . A -r 



B*B\ 



D 



+(-) 



-^B*B(A'A) if ^A'A > r_B«s, 
^ A* A(B* B) if < r_B*S- 



(66) 



(67) 



Let us analyze Eq. (|6J) for three particular cases. 

(i) If Ya-a >F b . s is valid, then from Eq. (||) it fol- 
lows that TA*Ato = ln4, and Eq. ( |64| ) reduces to 



Wl 



\K-b*a\ 2 /^a*a ■ 



(68) 



(ii) For Ta*a=^B'B the relation T B * B t Q = 2 — \/2 
holds. Thus, Eq. (p4l) reads as 



u>i ~ \K B . A \ 2 2 (V2 - i) e -( 2 -^)/r 



(69) 



(hi) In the case where T b * b ^>Ta-a is valid, one finds 
F B * B t =ln4, so that Eq. (p4[) reduces to 



Wl 



\K- B - a\ 2 /L ' b * b 



(70) 



A typical example of the temporal evolution of P B (t) 
for weak resonant dipole-dipole interaction is plotted in 
Fig. g for the case (ii). The rate regime with linear time 
dependence is seen to be established after some short 
time interval where a ^-dependence is observed. 



0.002 - 



(63) a? 



0.001 




FIG. 1. The probability of finding atom B in the 
upper state as a function of time for two atoms situated 
at diametrically opposite positions outside a dielectric mi- 
crosphere. The parameters have been chosen such that 
V a *a = Vb-b = 1.07r , Ta-b = 0.04r , and S a *b = 0.06r , 
To being the single-atom decay rate in free space. The solid 
curve is exact, in accordance with Eq. (p4|), and the dashed 
curve is Ps(t) — fs(io) + Wi(t — to). 

Another possible way for determining an energy trans- 
fer rate is to infer it directly from the equations of motion. 
From Eq. (pj) [or Eqs. @ and (§§)] we find that 



Q BB = -^ B ' B Q BB + (JC B 'AQA B + c.c), 
QAE = -\ (r^'A + Tb'b) QAE 

+ K,a*bQbb + K-*b*aQaa , 
where the notation 



Qaa 1 



= (U a \q\Ua>) = C A C* A 



(71) 



(72) 



(73) 



is used (u)a = wa')- 111 t ne rate regime qbb (= Pb) varies 
(approximately) linearly with time. This may happen for 
times where qab may be regarded as being (quasi-sta- 
tionary. By setting qab = in Eq. ([72]) and substituting 
the resulting expression for qab into Eq. (frD), we obtain 



Qbb — 



B'B 



2K,A'bK-B'A 
Taa» + ^b*b 



4|/C 



B*A\ 



Fa* a + Tb' 



c.c. 



■ QAA ■ 



QBB 



(74) 
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In the rate regime, the small first term in Eq. ( f74[ ) may be 
neglected, and in the second term qaa may be replaced 
with an appropriately chosen P A °^ = g AA ■ The energy 
transfer rate obtained in this way is then given by 



W2 



4|AC 



B*A\ 



>(0) 



B*B 



(75) 



The interest in energy transfer processes has stemmed 
from their key role in a wide variety of both biological and 
nonbiological systems, where the species participating in 
the energy transfer are typically (not necessarily identi- 
cal) molecules with dense manifolds of vibronic states. 
Commonly, a rate regime is considered and the trans- 
fer rate is calculated by means of Fermi's golden rule 
in second-order perturbation theory with regard to the 
molecule-field interaction, 



(76) 



where Wfi is proportional to the absolute square of the 
second-order interaction matrix clement and the energy 
conserving S- function S(ujf —u>i), and the sum runs over 
the continuum of initial (i) and/or final (/) states. In 
particular, short-distance energy transfer, where the in- 
termolecular coupling is essentially static, has been well 
known as Forster transfer pl[ . Later on energy transfer 
over arbitrary distances has been considered (see, e.g., 
Refs. HQ and references therein). The energy transfer 
rate in the presence of dispersing and absorbing material 
surroundings has been calculated in Ref. p4[ on the basis 
of a Hamiltonian of the type given in Eq. (13) . 

For two two-level atoms the single-transition probabil- 
ity per unit time uiif corresponds to Pb- Starting from 
the effective Hamiltonian (|2^) and evaluating Eq. ( |33| ) 
(for Cb) in first-order perturbation theory, we derive, on 
making the standard long-time assumption, 



P B = 2it\ICb*a\ 5(lu a ~uj b ) 



(77) 



(Appendix |^), which is fully consistent with the result, 
obtained in the standard second-order perturbation the- 
ory on the basis of the original (fundamental) Hamilto- 
nian (J2j) fyj . The only new feature is in the current 
treatment the medium induced atomic frequency shifts 
are taken into account. Note that the transition matrix 
element is determined by the full Green tensor. 

For establishing a rate regime, it is necessary that, ac- 
cording to Eq. (|76[), a continuum of initial and/or final 
states is involved in the transition. In the two-atom prob- 
lem at hand, the continua of states are obviously provided 
by the atomic level broadening due to the spontaneous 
decay. Let va and vb be the frequencies of the con- 
tinua associated with atom A and B respectively and 
£a(va) and £b(kb) the densities of the respective con- 
tinua. From Eqs. ( f76| ) and (77) the interatomic energy 
transfer rate is then expected to be 



ir = 27, I dv A J &vb£,a{va)£,b{vb) 
x Pa(va)\ICb*a\ 2 S{va 



vb) 



(78) 



Assuming Lorentzian line shapes 



1 



r A ,. A ,/2 



ir(v A > -Ca A >) 2 + (T A >*A>l2y 



(A' = A,B) 



(79) 



and making the assumption that p A (v A ) and |/Cb*a| 2 as 
functions of frequency slowly vary over the lines, we may 
perform the integrals in Eq. ((78]) to obtain 



w 



4|/C 



B*A\ 



Ta*a + r B » B 



Pa&a), 



(80) 



which exactly agrees with the rate wi given by Eq. J7q), 
if the assumption is made that 



p A {to A ) ^ P 



(o) 



(81) 



Finally, let us identify p A (ui A ) in E q. (p0| ) for the per- 
turbative rate w with P A (to) [Eq. (|5g) together with 
Eq. (|6^)] and compare the result with the nonperturba- 
tive rate Wi calculated from the exact temporal evolution 
of P B (t) [Eq. (g) together with Eq. @]. Noting that 
Pa (to) — exp(— TA'Ato), we find for the three cases con- 
sidered in Eqs. (p|) (f70|) the following results. 



(i) Ta'A^>^B*B- 



PA(to) 



w 



(82) 



(ii) Ta'A^B'B- 

P A (t ) ~ e-i 2 -^ , ^i~V2- 1-0.41. (83) 



(iii) Tb'B^^A'A- 

P A (t ) 



1, — 



w 



(84) 



Agreement between w\ and w is only observed in the first 
case, where Tb*b is sufficiently small. With increasing 
value of T b*b an increasing discrepancy between w\ and 
w is observed. 

The reason can be seen in the fact that there are two 
processes which simultaneously drive atom B: the en- 
ergy transfer from A to B and the spontaneous decay. 
Hence, Wi is actually the rate of both processes com- 
bined, not the "naked" rate of energy transfer. This ex- 
plains why w\ is typically smaller than w and why the 
discrepancy between them becomes more substantial for 
enhanced spontaneous decay of atom B. Note that the 
role of the spontaneous decay of atom B is twofold. On 







the one hand, it takes away the population of the upper 
state of atom B, thus diminishes w\ relative to w. On the 
other hand, it leads to an atomic level broadening which 
must carefully be taken into account for a proper eval- 
uation of the rate w according to Eq. (|7q). To roughly 
compensate for the first effect, the ratio wi/w may be 
multiplied by exp(TB*Bto). The values of the so cor- 
rected ratio are then ~ 1 for the cases (i) and (iii), and 
~ 0.74 for the case (ii). 



VI. POWER SPECTRUM 

Let us finally address the question of the influence of 
the resonant dipole-dipole interaction on the power spec- 
trum of light emitted by two identical atoms that initially 
share a single excitation. As is well known, the (physical) 
spectrum can be obtained by a Fourier transformation 
of the two-time correlation function of the electric-field 
strength in normal order (see, e.g., 25 1) 



S(r,wg,T) = 



die 



dti 



-iu>s(*2— ii) 



E(-)(M 2 )E(+)(rA) 



(85) 



where cus is the setting frequency of the (ideal) spectral 
apparatus, T is the operating-time interval of the detec- 
tor, and 



poo 

E (+) (r) = / du;E(r,w), 
Jo 



E ( - } (r 



E (+) (r) 



(86) 



(87) 



(Appendix |5]). Note that replacing Ci^A' — w) with 
■k8{u>a' —uS) would be too rough here. Equation (|88| ) re- 
veals that the emitted light is spectrally split into two 
asymmetric lines at uja T $a* b, with T± and \F a ± Fb\ 2 
being the widths and weights respectively. The line sep- 
aration is seen to be twice the dipole-dipole coupling pa- 
rameter. A system of two two-level atoms, one of them 
initially excited, is obviously equivalent to a three-level 
system with two upper dressed states |±). The doublet 
structure of the emitted-light spectrum can be easily un- 
derstood as a result of the transitions of the dressed states 
to the ground state. For the case of the atoms being in 
vacuum, similar results were found not only for the spon- 
taneous emission jl6| but also for the resonance fluores- 
cence (see, e.g., Ref. ^(| and references therein). 

As can be seen from Eq. (|S8|), an experimental obser- 
vation of doublet structure of the emitted light requires a 
delicate balancing act. The interatomic distance should 
not be too large to provide a reasonable level splitting, 
but it should not be too small to avoid \Fa\ = |Fs|, i.e., 
quenching of one of the two lines. It is worth noting 
that the presence of macroscopic bodies may facilitate 
the detection of the doublet, because it offers the possi- 
bility of realizing strong resonant dipole-dipole interac- 
tion even for interatomic distances much larger than the 
wavelength. 

Another interesting feature is that, according to 
Eq. (|50| ), either T + or can be much smaller than 
^a*a (= Ts's)- That is, the resonant dipole-dipole in- 
teraction can give rise to an ultranarrow spectral line, 
albeit each time at the expense of the other line of the 
pair. If, e.g., a single atom is placed sufficiently near a 
microsphere, its spontaneous decay may be suppressed, 
with the emission line being accordingly narrowed. Com- 
pared to the emission line of a single atom, one line of 
the doublet observed for two atoms being present may 
be further narrowed by several orders of magnitude l| . 



A. Weak atom-field coupling 



B. Strong atom-field coupling 



For weak atom-field coupling we derive, on basing on 
Eqs. (H), ©, ©, and (pi), 



S(r,u) S ,T -> °o) = ^ 



Au s + $A*B +tf+/2 



Aus - 5a*b +iY-/2 



where 



and {A' = A, B) 



Au! S = us-&a 



(89) 



3l : 

ire c- 



I dwImG(r,r A ',w)d^C(w A ,-w) (90) 
Jo 



For strong atom-field coupling, Eq. (^) changes to, on 
basing on Eqs. ([53|) and (|54]), 



5(r,w S) T^oo) 



(W A ±W B ) 



1 



Aw s ± 6 a *b + n±/2 + iAuj m /2 
1 

Auj s ±5 A -B- n±/2 + iAuj m /2 



+ 



i(FiTF B ) 



Au S T$A*B + iT T /2 



where (A' = A, B) 



(91) 
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w 



ImG(r,r A ',w m )dA' 



(92) 



(Appendix ^) . Here the upper (lower) signs again refer 
to the case where the state |+) (|— )) is strongly coupled 
to the medium-assisted field. Eq. ( pi|) reveals that due to 
the strong atom-field coupling, the doublet observed for 
weak atom-field coupling may become a triplet, with one 
of the lines of the doublet being split into two lines. These 
lines separated by Cl± have equal widths (which are solely 
determined by the width of the medium-assisted field res- 
onance) and equal weights. Note that their width and 
weight are different from those of the third line, which is 
closely related to a line of the doublet observed for weak 
atom-field coupling. From Eq. ( |9l|) it is also seen that, 
depending on the point of observation, this line or the 
strong-coupling-assisted doublet can be suppressed due 
to the interference effects. 



VII. SUMMARY AND CONCLUDING REMARKS 

We have studied the interaction of two-level atoms with 
the electromagnetic field in the presence of dispersing and 
absorbing material surroundings described by a spatially 
varying permittivity that is a complex function of fre- 
quency. Starting from the exact multipolar Hamiltonian 
(in electric-dipole approximation), we have derived re- 
duced density-matrix equations of motion that describe 
both the on-resonant interaction of the atoms with the 
medium-assisted electromagnetic field and the resonant 
dipole-dipole interaction of the atoms with each other 
via off-resonant atom-field interaction. The equations 
of motion apply to atoms in the presence of arbitrarily 
configured, dispersing and absorbing media and are not 
restricted to the rotating-wave approximation. Whereas 
the resonant dipole-dipole interaction is essentially con- 
trolled by the real part of the medium-assisted Green 
tensor, the relevant quantity for the resonant atom-field 
interaction is the imaginary part of the Green tensor. 
Both together determine the mutual interaction of the 
atoms. 

We have shown that when the differences between the 
atomic transition frequencies are small compared to the 
frequency scale of variation of the medium-assisted Green 
tensor, then the equations of motion can be derived from 
an effective Hamiltonian. All coupling parameters are 
again expressed in terms of the medium-assisted Green 
tensor. The use of the effective Hamiltonian may sub- 
stantially simplify the calculations in cases where iden- 
tical atoms are considered. Applying the theory to the 
two-atom case by assuming that initially the medium- 
assisted electromagnetic field is in the ground state and 
one of the atoms is in the upper state, we have considered 
the resonant energy transfer between the atoms and the 
spectrum of the emitted radiation. 

There are two energy-transfer channels in general: one 
channel through resonant dipole-dipole interaction medi- 



ated by virtual-photon creation and destruction and the 
other one through emission and absorption of real pho- 
tons. In particular for strong atom-field interaction, the 
(over a period averaged) energy transfer can be inhib- 
ited due to destructive interference of the two available 
transfer channels. When the dipole-dipole interaction is 
weak, the energy transfer from the donor atom to the 
acceptor atom becomes irreversible, and a transfer rate 
can be inferred from the (exact) temporal evolution of 
the excited-state population of the acceptor atom. Such 
a rate regime is commonly described by a transition prob- 
ability per unit time, such as the Forster transfer rate, 
which is calculated by means of Fermi's golden rule in 
second order perturbation theory with respect to the 
fundamental Hamiltonian. We have shown that the two 
methods essentially lead to the same expression for the 
transfer rate. In any case it is proportional to the abso- 
lute square of the medium-assisted (full) Green tensor. 
In contrast to molecules, where the necessary continuum 
of initial/final states is typically given by the vibronic 
states, in the case of atoms this continuum is essentially 
built up during the process of spontaneous decay. 

The results show that both the resonant energy trans- 
fer and the doublet spectrum of the emitted light ob- 
served for weak atom-field coupling or the triplet spec- 
trum observed for strong atom-field coupling can be con- 
trolled by the presence of macroscopic bodies. Clearly, 
the present analysis has left a number of open questions, 
on which future work will concentrate. In particular, the 
problem of energy exchange between atoms whose transi- 
tion frequencies must be regarded as being different with 
regard to the variation of the medium-assisted Green ten- 
sor, needs special emphasis. 
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APPENDIX A: DERIVATION OF 
Eqs. (§) and (^) 

In order to perform the second term in Eq. (|5|), 

F"(t) = -J2j o [[6{t)M(t)]MrA,w,t) 

+ ±\v A ,uj,t) [d{t),d A {t)]} , (Al) 
we formally integrate Eq. (^) to obtain 



u; 2 si(r,w) 



hiteo 
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x V / dt'd A (t')G(r A) r,uj)t 

A JO 



j(t-t') 



(A2) 



where ff rGC (r, u),t) evolves freely. Inserting Eq. (A2) into 
Eq. (Q), we derive 

E(r,w,i)=E free (r,a;,t) 

+ — ^ E / dt'e"^ C-'VGfr.rA.^dxff), (A3) 

7T£ C 2 ^ Jo 

where E frco (r, a;, t) is denned according to Eq. (Q) with 
ff re e(r, ui, i) in place of f (r, oj, i). 

Let us first restrict our attention to the single-atom 
case. Introducing slowly varying atomic operator 



we may rewrite Eq. (0) as 



d A (t) = d A a A (t) 



, — iuAt 



it 



d A a A (t) 



Jul At 



(A4) 



(A5) 



Inserting Eqs. (|A3|) and (|A5|) into Eq. (Al), we obtain 



i ft {."ac 2 r r 



■ d^(t')e _l( " +WA)(t_ *' ) e iWAt 



- Im G(r A , r A , w) [d A a j4 (t')e 4(a,+ ^ )(t_ *' ) e-' i "- 4 * 

+ d^^(t')e i( "^ A)(t -*' ) e^- 4t ] [6(i),d A (*)]}, (A6) 
where 

-^AfreeC*) = ' 



dc^ {[6(<),dx(t)]E free (r x ,w,t) 
+ EL(r A ,w,t) [0(t),dA(*)]}. (A7) 



Confining ourselves to resolving times large compared 
with l/(u> — uA we may regard the exponential func- 
tions in Eq. ( |A6| ) as being rapidly varying with time. 
Recall that range of integration only covers off-resonant 
frequencies. In the time integrals, we may then replace 
the slowly varying operator a A (t') by a A (t) (Markov ap- 
proximation). In the spirit of the coarse-grained time 
averaging mentioned, we may further make the replace- 
ment 



(A8) 



[((x) = ttS(x) + iV/x] and drop the terms where a A or 

a A appears twice. In this way, Eq. ( |A6| ) approximately 
reads as 



1 f"°° UJ 2 

/ duj — d* A lmG(r A ,r A ,uj)d A 

hire Jo c 2 

x {c[-(w + w 4 )j [6(t),a A (t)\ a\{t) 
+ ((uj a -lu) [d(t),4(t)j a A (t) 
- ((uj + uj a ) a A {t) \6(t),a\(t) 



C(io-Lo A )a A (t)[d(t),a A (t)\y (A9) 



Substituting Eq. (A9) into Eq. (g) and taking into ac- 
count that, because of u>± ui A y^Q, the (5-function parts 
of the (^-functions do not contribute to the frequency in- 
tegrals, we make the approximation 

C(w ± u A ) -> iV/(u ± uja) (A10) 
and arrive at 



o = — 



0,H 



AS 



TP" 

r A free • 



(All) 



Here, Has is the single-atom system Hamiltonian de- 
fined according to Eq. with the "naked" transition 
frequency uja being replaced by the shifted frequency 
uja — uj a — $A*Ai where the frequency shift 5a* A is given 
by Eq. (|ll|) together with Eq. @. 

Let us turn to the multi-atom case. Substituting the 
formal solution (A3) into Eq. (Al), we now obtain a 
double-sum over A and A' , 



(A12) 



A, A' 



The terms F A (t) are treated as described above to give 
the shifted transition frequency for each atom. In the 
terms F AA , (t) (A ^ A') we then take into account the 
level shifts and introduce the slowly varying atomic op- 
erators according to 



d A {t)=a A {t)e iQAt . 



(A13) 



By repeating for the terms F AA ,(t) the same procedure 
as in the single-a tom case, it is not difficult to prove that 
[in place of Eq. ( A9)] the result is 



1 



ui ( 

dw— jd A ImG(r A ,r A ,,uj)d A , 



xC[-(w+5U0] [oW.^wj^W 

- d* A Im G(r A , r A , , uj) d A >C(Q A > -«) [d(t), a A (t)\ a A , (t) 
d* A Im G(r A , t a ,,u) d A ,((uj+uj A ,)a A ,(t) [6(i), a A (t) 

- d A Im G(r A , t a ,,u) d* A ,((uj-uj A ,)a A , (t) [o(t), cta (*)]}. 

(A14) 



12 



In Eq. (A14) we now make the approximation 

C(v±WA>)^iV/(w±&A'), (A15) 



[cf. Eq. flAlCD l and insert the resulting expression for 
Pj A , it) I into Eq. ( |A12[ ). We eventually combine Eqs. (§) 
and ( A12| ) and obtain, on again dropping off- resonant 
terms, 



= -nl°> H * 



+ i 



A' A' 



A. A' 



+ S AA'* a A' 



TP" 
^fre 



6,4 



0,(7,4 



0,a A 



6,4 



where 



x free / y ± A free ' 



}, (A16) 



(A17) 



and Saa' i s g iven by Eq. ( fl3| ) together with Eq. (14). 
From Eq. (A7) it is seen that when the off-resonant free 
field is in the vacuum state, then 



(^r'ee)=0 



(A18) 



is val id, an d the expectation value of the operator equa- 
tion ( |A16| ) just yields Eq. (|). 

To prove Eq. ([n]) , we use the cyclic properties of the 
trace and rewrite Eq. (H) as 



dt 



(6) = Tr 







H 


H s ,g 



+ 1 



[ 6 A'A> {f A v A 'Q-o A 'Qo\ 

A,A> 

+ 5 AA'* {°Ao\,Q- b\,QCTA 
+ S AA'* (°AQVa> - Q^A'^A 

+ S A*A> (Aa&A' ~ QPA'Oa 



6 . (A19) 



The relationship ([l8|) obviously implies the relationship 



Tr 



dg{t) 
dt 



0(0) 



= Tr 



£(0) 



dO(t) 
dt 



(A20) 



Since O is an arbitrary system operator, from Eqs. ( A19) 
and (A20) it then follows that Eq. (|l9|) must be valid. 



APPENDIX B: DERIVATION OF Eq. fl20| ) 

For weak atom-field coupling, the on-resonant term 



= JE dLU i \0(t),d A (t)] E(r A , w, t) 
n A Jo lL J 

+ ±\r A ,U,t) [6(i),d A (*)]}, (Bl) 



in the Hamiltonian H$ in Eq. ( A16| ) can also be treated 
within the approximation scheme outline d in Appendix 
[A| for the off-resonant term F"(t) [Eq. (Al)]. Writing 
F'(t) as 



F'{t) 



Ffre 



A.A 



F A A<(t) 



(B2) 



where F/ ree is defined according to Eqs. (A7) and ( |A17| ), 

but with Jp 00 dcu ... in place of J Q °° dcu . . ., and repeating 
the steps leading to Eq. ( A14[ ), we obviously arrive at an 
equation for F' AA , (t) which again looks like Eq. flAT^ ), 
but again with J °° dui ... in place of J Q °° dcu .... Note 

that the equation for F AA ,(t) also applies to the case 
A = A'. Since the frequency integrals now run over the 
resonance region, only the 5-function parts of the £- 
functions do effectively contribute to the frequency in- 
tegrals, i.e., we may approximate the £-functions by S- 
functions, 



((u> ± Cja') — ► ttS(ui ± Cja')- 



(B3) 



Inserting Eq. (B2) [together with the equation for F' AA , (t) 

and Eq. ( p33j )] into the expressi on fo r Hg in the first term 
on the right-hand side of Eq. (A16), we derive 



O 



a 

A.A' y 



A, A' 



A* A' 



- 5 AA>* & 



0,(TAz 



0,a\ 



6,4 



0,a A 



F 



free 



o\4' + r AA >*a A , 



0,a A 



*a.+5 + aa , 



+ 5\, A ,a A > 



0,a A 



6,4 



4' 



where Taa 1 is given by Eq. (21), and 

-Ffrec = -Ffrec + ^free" 



(B4) 



(B5) 



Assuming that the medium-assisted electromagnetic free 
field is in the vacuum state, thus 



(B6) 



(F free ) = 0, 



and recal ling the procedure that has led from Eq. ( |Alq ) 
via ( |A19D to Eq. ([l9|), we see that Eq. (|B4|) just leads to 

Eq. W 
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APPENDIX C: DERIVATION OF Eq. (|77|) 



To calculate P B (t) = \C B (t)\ 2 in lowest-order of per- 
turbation theory with regard to Eq. (|3^) , we replace C A 
and Cb on the right-hand side of the equation for C B by 
their respective initial values 1 and 0, and integrate both 
sides of the equation with respect to time to obtain 

C B (t)=i5 B * A I <tf e i @ a - aA) * 



+ dt' dt" / duK B , A (t',t";u)). (CI) 
Jo Jo Jo 

Using Eq. together with the relationship 



At'-t")e{t' -t") 



1 

'2ttI 



-iu'{t'—t") 



do/- 

oo Uj' -UJ + tO + 



(C2) 



we may rewrite Eq. (CI) as 
C B (t)=iS B , A 2n5^(u> A -uj B ) 



+ 



the. 



do/ 



J- 



,2 A* 



duj ■ 



d* B JmG(r B ,r A ,uj)d A 



uj' — uj + i0 + 



xS w (uj' -u) B )5^(u) A -u'), 



where 



<5 (t) M = ^ fdt'e 
27r Jo 



ILOt 



(C3) 



(C4) 



The function A (a/) = 6® (uj' - u> B ) 6® (Q A - uj') in Eq. 
(C3) is essentially different from zero in an interval 
around uj a ~uj b , the extension of which is of the order 
of magnitude of 1/t. For sufficiently long times t, it is 
reasonable to assume that the function 



R(uj') = 



uj 2 d* B ImG(r B ,r A ,uj)d / 

c 2 uj' — uj + i0 + 



(C5) 



is slowly varying on the frequency scale of variation of 
A(u/). The integral over uj' can then be performed sep- 
arately to yield 



duo' <5 (t) (J - Cj b ) 5 {t) (Cj A - uj') = 5 {t) (Cj A — u> B ), 

(C6) 



and Eq. (C3) takes the form of (t— »oo) 
2 



C B (t) = 



2m8 



B' A 



iheo 



uj' 2 d* B Im G(r Bl r A ,uj)d A 

c 2 uj a — uj + i0+ 



6^(uj a -uj b ). (C7) 



To evaluate the frequency integral in Eq. (|C7|), we com- 
plete it by adding the corresponding off-resonant part 
and subsequently subtract it. The integral J °° duj . . . 
can then be approximately evaluated by extending the 
lower limit to — oo and using contour integral technique 
in a similar way as in Ref. jlj]. The remaining inte- 
gral J °° duj . . . can again be treated as a principal- value 
integral to give —2iriS B * A [cf. Eq. (0)], which (apart 
from the quantum correction) cancels out the first term in 
Eq. flC7|). Taking into account that for t — >oo, \S (t \uj)\ 2 
behaves like 



\ 5 W(uj)\ 2 = ±-S(uj), 



(C8) 



we eventually arrive at 



P B (t) = \C B (t)\ 2 = 



~,2 



2tt 



h 2 \e c 2 



x \d B G(r B ,r A ,cu A )d A \ 2 S(uj A - Q B )t, 
which together with Eq. (37) leads to Eq. (|77|). 



(C9) 



APPENDIX D: DERIVATION OF 
Eqs. dH) and (|l|) 



Let us briefly outline the derivation of Eqs. ( pq ) and 
( |9l| ) (for more details, see [[Lsl). In the Schrodinger pic- 
ture, the two-time correlation function in Eq. ( |S5| ) can be 
given by 



E(-)(r,t 2 )EW(r,t 1 ) / 
= U(t 2 ) E<-'(r)e-*" (i! -' l)/t E( + '(r) (Dl) 



Calculating it by using Eqs. (|30|), (]3J), (J86J), (|87|), and 
(Q), and inserting the resulting expression into Eq. 
we derive 



S(r,oj s ,T) 



duj ■ 



V f dti ^-^/WcaM 
A Jo L Jo 



7reo<r 



ImG(r,r^,w)d4e 



U>A)(tl— t') 



(D2) 



For the two-atom system under consideration, in 
the weak-coupling regime, C A ( B )(t') may be replaced 
with CA( B )(t\), where CU(B)(ii) is given according to 
Eqs. © and ©, 

C A (B)(t) = \ 

X [ e (- r +/2+«A*B)t + (_v(-r_/2-i5 A « B )tj _ / n3 N 
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(Ta*a — ^b*b and Ka*b — K>b*a)- The ^'-integral may 
then be regarded, in the long-time limit, as being 
C,{u>A' ~ UJ )i so that the factor F^' [Eq. (|90|)] can be put 
in front of the ii-integral. Now the ti-integral can be 
performed to obtain Eq. ( §81 ) (T— *oo). 

For strong a tom -field coupling, we use Eq. ([l5|) and 
express in Eq. (D2) CA(B){t) m terms of C±(t), 



S(r,wg,T) = - 



1 ' 

/ di ie 4( 
Jo 



U$—u>A)ti 



tl 



dt'C + (i')e i<W *' 



dw ■ 



JO Te C 

x pmG(r, r j4 ,cj)d j4 +ImG(r,r B ,w)d B ] 

— WA)(ii-*') 



x e 



+ / dt'C_(t')e" 1 ^ 



i<5,,» H t' 



DC. 



dw 



o Jo teoc* 

x [Im G(r, , di - Im G(r, r B , w) d B ] 



-i(a>— ou)(*i — *') 



(D4) 



Recall that either the state |+) or the state |— ) is strongly 
coupled to the medium-assisted electromagnetic field, but 
not both at the same time. The state which is weakly 
coupled to the field can be treated in the same way as 
above. For the strongly coupled state, we again assume 
a Lorentzian shape for the field res onanc e (w m , central 
frequency; Au m , width; cf. Section IV B ) and evaluate 
the w-integral. Taking C±(t') from Eq. (|53|), we then can 
evaluate the remaining time integrals to obtain Eq. j9l| ) 
together with Eq. @ (T->oo). 
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